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Zusammenfassung 



Neutrinophysik hat in diesem Jahrzehnt sehr an Popularitat gewonnen, nicht 
zuletzt aufgrund der Messung von Neutrinooszillationen. Allerdings lassen 
sich nicht gleichzeitig alle Oszillationsexperimente vollstandig durch das mit 
Neutrinomassen erweiterte Standardmodell erklaren. Eine mogliche Losung 
bietet die Einfuhrung eines vierten Neutrinos, das jedoch aufgrund der Z° 
Resonanzbreite nicht schwach wechselwirken darf. 

Ein solches steriles Neutrino konnte aufgrund seiner gravitativen Wech- 
selwirkung sowie der Neutrinooszillationen einen EinfluB auf die primordiale 
Nukleosynthese haben. Dieses Phanomen ist schon in einer Vielzahl von 
Veroffentlichungen behandelt worden. Es hat sich jedoch als sehr schwierig 
herausgestellt, das nicht-lineare Differentialgleichungssystem, dafi die Neu- 
trinooszillationen im friihen Universum beschreibt, zu losen. 

Wahrend bisherige Publikationen sich meistens auf numerische Rechnun- 
gen beschrankten, untersuchen wir das System analytisch, wobei wir als 
einzige Naherung die Impulsverteilung vernachlassigen. Dabei achten wir 
besonders auf das Verhalten des Systems kurz nachdem die Resonanztem- 
peratur unterschritten worden ist. Im Einklang mit anderen Publikationen 
sehen wir einen exponentiellen Anstieg in der Neutrinoasymmetrie. Wir 
beweisen, daB fur einen signifikanten Bereich der Mischungsparameter 5m 2 
und sin 2 29 diese Neutrinoasymmetrie zu oszillieren beginnt und sogar das 
Vorzeichen wechselt. Damit werden numerische Kalkulationen bestatigt, die 
diesen Effekt schon fruher entdeckt haben, deren numerische Stabilitat je- 
doch angezweifelt wurde. 
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Chapter 1 
Introduction 



Neutrino physics has become very popular during the past decade. The main 
impetus in this field has come from the strong improvements in detecting 
neutrinos; there is now strong evidence from two natural neutrino sources, 
the solar core |]] and the Earth's atmosphere @, that neutrinos of a certain 
flavor disappear on their way to the detector. Furthermore, a neutrino beam 
produced in the laboratory has been found to change flavor in the LSND 
experiment @. 

The most natural explanation for these anomalies is achieved by extend- 
ing the Standard Model with neutrino masses. Together with the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix for the lepton sector, the masses imply 
flavor mixing; neutrino oscillations have been born |4]]. Such oscillations 
can be fully parameterized by the differences between the squared masses 
5mfj = m 2 j — mf and the mixing angles 9ij. Also, a possible phase in the 
CKM matrix could induce CP violation. With these neutrino oscillations, 



the experiments can be explained by the solutions presented in Table [TTT . 

However, with this straightforward extension of the Standard Model it is 
not possible to explain all three experiments simultaneously, simply because 
the trivial condition 

J2 Sm ij = ( m 3 - m l) + ( m l ~ m i) + ( m i - m l) = (1-1) 

cannot be fulfilled by the data. Hence, many interesting new models have 
been suggested to explain the three different experiments. One of the most 
appealing possibilities is the existence of a fourth neutrino. Since the Z° 
decay width || excludes more than three standard neutrinos, this new flavor 
must be inert with respect to the standard weak interaction. 
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Experiment Favored Channel \5m 2 \ [eV 2 ] sin 2 26 



Solar 

Vacuum u e — > anything (0.5-8) x lCT 10 0.5-1 

MSW (small angle) u e -> anything (0.4-1) x 10~ 5 lO^-lO" 

MSW (large angle) u e -> ^ or z/ T (3-30) x 10" 5 0.6-1 

Atmospheric — > z/ r (1~8) x 10~ 3 0.85-1 

^ -> i/ 8 (2-7) x 10~ 3 0.85-1 

LSND -» z/ c 0.2-10 (0.2-3) x 10 
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Table 1.1: Results from neutrino oscillation experiments j^j. The values are 
nominal 2a ranges. 

Although the hypothesis of such a sterile neutrino may seem highly spec- 
ulative, its possible existence is the most far-reaching implication of the 
current experimental situation. It is certainly worthwhile to investigate its 
consequences! Sterile neutrinos would have a strong impact on certain as- 
trophysical environments due to their mixing with active neutrinos. These 
effects would become especially important under extreme conditions, like in 
core-collapse supernovae |7| and in the Early Universe. 

It is the impact of sterile neutrinos in the Early Universe, and in particular 
the outcome of the primordial nucleosynthesis, that has been the subject of 
an intense recent debate. By incorporating sterile neutrinos into the well 
understood mechanism of standard Big Bang nucleosynthesis (BBN), two 
frontiers of physics would be connected in one phenomenon, making this a 
powerful tool for predictions. 

The production of primordial nuclei depends sensitively on the expansion 
rate of the Universe, which in turn depends on the number density of sterile 
neutrinos. This simple correlation is well known ||; the actual challenge 
is to deduce how much the sterile neutrino sector is populated. Since the 
only connection between sterile neutrinos and ordinary matter, apart from 
gravity, is given by oscillations, this phenomenon must be the key for creating 
a relation between sterile neutrino parameters and Big Bang nucleosynthesis 
predictions. 

In the beginning of the 1990 's, the number density for sterile neutrinos 
at the time of primordial nucleosynthesis was calculated depending on the 
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parameters for oscillations between an active and a sterile neutrino, 5m 
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and 9 || [T(J. The idea was to constrain the allowed parameter space, since 
strong observational results (which have weakened since) set an upper bound 
on the number density of the v s . In these calculations, the influence of the 
thermal plasma on the mixing parameters was included, but the small CP 
asymmetric contributions were neglected. The constraints were so strong 
that they clearly excluded the <-> u s solution for the atmospheric neutrino 
anomaly. 

In 1995, Foot, Thomson and Volkas JTTJ] found an interesting effect by 
including the CP asymmetric contributions; the asymmetry could induce dif- 
ferent population rates for the sterile neutrinos and antineutrinos, or equiva- 
lently, different depopulation rates for the active neutrinos and antineutrinos. 
This would change the neutrino asymmetry and thus would have an influence 
on the population rates. They found that this back-reaction could amplify 
an initially small CP asymmetry by several orders of magnitude for a large 
range of parameter space. (Actually, this back-reaction had been discussed 
much earlier by Barbieri and Dolgov ||10|| , but they erroneously found the 
effect to be small.) 

The implications from this new effect were manifold. The most interesting 
conclusion was given by a model in which v T <-> u s oscillations created a large 
v T asymmetry, which in turn suppressed <-> u s oscillations [12-17]. In such 
a model, the atmospheric neutrino anomaly could be explained by <-> u a 
oscillations without the sterile neutrinos coming into thermal equilibrium. 
Another thoroughly discussed implication treated the impact of a large u e 
asymmetry on primordial nucleosynthesis |13, pl| the large chemical 



potential /x e would change the neutron-to-proton ratio and thus the nuclear 
abundances. 

Clearly, sterile neutrinos can only be included correctly into BBN if the 
mechanism producing large neutrino asymmetries is revealed. Unfortunately, 
this has proven to be a very difficult task. A series of papers have been pub- 
lished on this subject by several groups [12-36], but the situation remains 
unclear. The main problem is that the system of differential equations de- 
scribing the mixing is very complex: one has to treat neutrinos of different 
momentum separately, and one especially has to take care of the non-linear 
terms in the equations. Therefore, most works tried to solve the problem 
with numerical calculations. Thereby, they often used either the adiabatic 
approximation, which simplifies the equations for a given momentum, or they 
neglected the momentum distribution. Only during the past year some au- 
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thors claim to have solved the full system of differential equations numerically 



14, 15, 20 



Due to these difficulties in solving the system, the different works present 
contradictory solutions: when neglecting the momentum distribution, the 



neutrino asymmetry shows an oscillating behavior 21, 22[ . The system then 



ends up with a calculable value of the v asymmetry, but with unpredictable 
sign. Such a scenario would introduce domains in the Early Universe with 
different signs of neutrino asymmetry | |23|| . On the other hand, applying the 
adiabatic approximation [|13|] is questionable since the neutrino oscillations 
are not adiabatic close to the resonance for the parameter space of inter- 
est. Here the neutrino asymmetry does not oscillate. Since the numerical 
calculation using the full system of differential equations is very CPU-time 
consuming, one also has to question whether the calculations are done with 
sufficient accuracy and whether the numerics are stable at all. After all, 
non-linear systems tend to show chaotic behavior. Besides, different results 
achieved by these calculations are contradictory on the point of oscillating 
neutrino asymmetry [|15|, |23|j . 

All these numerical works have in common that they predict a similar final 
absolute neutrino asymmetry (if effects of neutrino domains are ignored) of 
orders 10 _2 -1. But even this outcome has been questioned very recently in 
an analytical work by Dolgov et al. [2~3|, who end up with a final asymmetry 
several orders of magnitude lower than the former results. Thus, the only 
point on which all works agree is that oscillations between sterile and active 
neutrinos in the Early Universe have the potential of creating a neutrino 
asymmetry at least of order 10 -5 , a number which is still large compared to 
the baryon asymmetry, rj = O(10~ 9 ). 

The main problem is to get a grip on the exact behavior of the non-linear 
system. Our work will try to bring some clarity into this subject. Most of the 
former papers were strongly based on numerical calculations. In this paper, 
the discussion will be based on a thorough analytical treatment, which will be 
supported by simple numerical calculations. We will describe the evolution of 
the system, neglecting the momentum distribution, and compare our results 
with a numerical solution for the parameters (5m 2 , sin26* ) = (—1 eV 2 , 10~ 4 ). 
It should be the job of future investigations to include the effects from the 
momentum distribution. 

In Chapter 2, those aspects of the Early Universe and Big Bang nucle- 
osynthesis will be summarized which are important for our system. Chapter 
3 contains all relevant information about neutrino oscillations in matter. In 
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Chapter 4, we analyze the system in the simple two flavor case. After intro- 
ducing initial conditions, we prove that our system creates oscillations of the 
lepton asymmetry. In Chapter 5, we summarize our findings. 
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Chapter 2 

Big Bang Nucleosynthesis 



We study the influence of neutrinos, in particular their effective number of 
flavors and the u e chemical potential, on Big Bang nucleosynthesis (BBN). 
To this end we first describe the expansion of the Early Universe before and 
during BBN. We define number density and lepton asymmetry. Next, we 
give a short summary of BBN, concentrating on the main product, 4 He, and 
we estimate the influence the neutrino sector can have on its abundance. 
Finally, we summarize the observational results. 

2.1 Dynamics in the Radiation Epoch 

In order to discuss the influence of neutrinos on BBN we need to introduce 
some general concepts pertaining to the relevant cosmological epoch. Natu- 
rally, we only need to look at the Universe at times before and during BBN, 
i.e. at temperatures above 0.1 MeV. On the other hand, for the mixing pa- 
rameters we will consider, neutrino oscillations are strongly damped at tem- 
peratures much higher than the QCD phase transition scale Aq C d of around 
200 MeV. The temperatures of interest are therefore between 200 MeV and 
0.1 MeV, corresponding to cosmological time scales from 10~ 5 s to 3 min. 

The Universe during this epoch is flat, homogeneous and isotropic, and 
can therefore be described by the Robertson- Walker metric. The Friedmann 
equation for these conditions is simply || 

H\t) = ^p(t), (2.1) 

where H(t) = R(t)/R(t) is the Hubble parameter or expansion rate, R the 
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cosmic scale factor, G = m~ L 2 Newton's constant, and m p \ = 1.22 x 10 19 GeV 
the Planck mass. Because the energy density p is dominated by radiation, it 
scales as i?~ 4 , which upon integrating ( |2.1| ) gives H = \ t~ 1 . 
The radiation density is usually expressed in the form 

P = ^9*T\ (2.2) 

where T is the photon temperature and g* the total effective number of 
relativistic degrees of freedom, 

<?*= E 9,(j) 4 +l E 9i(^)- (2.3) 

i=bosons i=fermions 

Here, gi is the number of internal degrees of freedom of particle species i and 
Tj is its temperature. If we insert (|2.2|) into ( |2.1|) and use t = 1/2H, we 



finally derive a relationship between time and temperature: 




"° S -V»^ = 0.301,^^. (2-4) 



3 ^* rp2 ^* J^2 



As another consequence, we note that T oc R~ l as long as g* = const. 

In the range 1 MeV <^ T <> 100 MeV, the only particles which contribute 
significantly to the radiation density are photons (g 1 = 2), electrons and 
positrons (g e = 4), and three left-handed neutrino families with g v — 2 each. 
For all of them Tj = T applies as long as they remain in thermal equilibrium 
so that g* = 10.75. Any additional neutrino species would add another 7/8 
per internal degree of freedom. For T J> 100 MeV, g* is higher due to the 
presence of muons and pions, and for T ^ Aqcd ~ 200 MeV many gluon 
and quark degrees of freedom are excited. For T <^ 1 MeV, the electrons and 
positrons become non-relativistic and annihilate, reducing g* by (7/8)g e = 
7/2. This effect heats the photons (the neutrinos are already decoupled) so 
that T v /T is reduced to (4/11) 1 / 3 , leading to g* = 3.36. 

In our subsequent discussion we will ignore these effects and always use 
the value g* = 10.75, implying 

1 T 2 

H=-t- 1 = 5.5 . (2.5) 

2 m p i 

This is a good approximation as long as the oscillation parameters are such 
that the crucial events happen in the range 1 MeV <^ T ^ 100 MeV. 
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2.2 Number Densities and 
Fermion Asymmetries 

Because we want to study the influence of neutrino number densities and 
matter-antimatter asymmetries, we introduce appropriate measures for them. 
It is convenient to normalize these quantities to the photon number density 



C: 



n 7 = 2 4T 3 , (2.6) 

where ( is the Riemann zeta function with £3 ps 1.202. 

Fermions may have non- vanishing chemical potentials /ij which affect their 
number densities. In addition, a non-vanishing /ij implies a CP asymmetry 
for species % which we parameterize as 

Ai = (2.7) 

where refers to the antiparticle density. 

At temperatures below the QCD phase transition virtually no hadronic 
antimatter exists. Since there are no baryon-number violating interactions at 
this late epoch, the total baryon number is conserved. Therefore, the number 
density of baryons, n B , scales as -R" 3 or T 3 . The baryon asymmetry 

A B = nB ~ % « ^ (2.8) 

is thus constant under the cosmic expansion. 

The commonly used present-time baryon asymmetry r\ is related to our 
A B through 

7] = -^A B = 1 to 6 x 10~ 10 . (2.9) 

The difference is due to photon heating when electrons and positrons annihi- 
late, an effect which modifies all asymmetries. We will always take A { to refer 
to the fermion asymmetry of species % at the epoch before e + e~ annihilation. 

For T < Aqcd, almost all baryons are either neutrons n or protons p, 
which implies that n B = n n + n p . The non-trivial evolution of n n and n p will 
be discussed in the next section. 

We will assume that the Universe is charge neutral. Neglecting muons 
and pions, we therefore use 

A C = A P (2.10) 
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for the electron asymmetry. 

Turning to neutrinos, we can safely neglect their mass. Therefore, their 
number densities are given by 




(2.11) 



where E = p = |p|, f(p, T, jj,) = [e( p ~^/ T + l]" 1 is the Fermi-Dirac phase- 
space distribution function, and fi is the chemical potential. For very small 
asymmetries (/i <C T) we find 

n v + n v = -n 7 + C((/i/T) 2 ) (2.12) 

and 

K = | + O (WT) 3 ) » 0.68 £ + O ((/i/T) 3 ) , (2.13) 

where we have used that the chemical potential for antineutrinos is p, — —[/,. 
Of course, these equations are only valid as long as the neutrinos are in 
thermal equilibrium. 



2.3 Helium Abundance 

Big Bang nucleosynthesis is a rather complex system depending on a number 
of parameters, including the baryon asymmetry rj and the effective number of 
neutrino families iV^ ff , or more generally, g*. This system has been analyzed 
very thoroughly, including exhaustive numerical calculations. Here we will 
give a short summary, concentrating on the production of 4 He. 

The outcome of BBN will be abundances of the different species (A, Z) 
with Z protons and A — Z neutrons. We will express these abundances as 
the mass fraction contributed by the species, i.e. 

X„ = 42±» (2.14) 
n B 

where ua,z and ne are the number densities of the nuclear species (A, Z) and 
all baryons, respectively. 

In strict thermodynamic equilibrium, the abundances are given by nuclear 
statistical equilibrium (NSE) so that 

(X A , Z ) NSE oc V A ~ l exp (-^) , (2.15) 



11 



where B^.z is the binding energy of the nuclear species (A, Z). We see that 
the NSE abundances are strongly suppressed by the small baryon asymmetry 
1] m 10 -9 . This suppression is compensated by the exponential factor in 



27T5|) at low temperatures T < B A ,z- For 4 He, (X He )NSE ~ O(l) at T » 
0.3 MeV. Heavier nuclei have significant NSE abundances only at even lower 
temperatures. 

The NSE abundance for an element X will freeze out when the rates 
T(ab <-> X) of producing it from the lighter elements a and b become smaller 
than the expansion rate H. We have 



T(ab — > X) oc n a nb exp 



k \ 1/3 ' 



T 



MeV 



(2.16) 



where fc « ZfZ^A a A b /(A a + A b ) and T Me v = T/ MeV. The first two factors 
are the number densities of the nuclear species a and b, respectively, the 
last factor represents the Coulomb-barrier suppression, which increases with 
Ai and Z^. The freeze-out temperatures increase with the Coulomb-barrier 
suppression. 

As the temperature decreases, the NSE abundances increase, but at the 
same time the nuclear reactions begin to freeze out. Therefore, heavy nuclei 
are not produced during BBN because they freeze out long before their NSE 
abundances have become significant. 

Apart from traces of other nuclei, BBN produces primarily 4 He so that it 
is a good approximation to assume that all neutrons end up in 4 He. Then 

v _ Y ~ 4nHe 4 K/ 2 ) 2(n D /n p ) BB N (0 1 7 ^ 

r = X Ue ~ = ; = — — -, j — r • (2.17) 

n B n n + n p 1 + (n n /n p )BBN 

Therefore, the all-important helium mass fraction Y depends primarily on 
the n/p ratio at the time when 4 He freezes out of NSE, which happens at 
^bbn = l - 3 minutes. 

To find (n n /n p )BBN, we have to follow the evolution of the n/p ratio from 
the beginning. At temperatures above 1 MeV, reactions of the type 

n + z/ e < — >p + e, (2.18) 

maintain chemical equilibrium with the rate T oc GpT 5 . Therefore 




— — - — 



(2.19) 
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where Q = 1.293 MeV is the mass difference between neutrons and protons, 
and fj, e and \i Va are the chemical potentials of the electrons and electron neu- 
trinos, respectively. Surely, \i e can be neglected since fi e /T « r\ ~ 10" 9 . For 
the moment we will set the neutrino asymmetry to zero which is equivalent 
to /i„ e = 0. Then the n/p ratio depends only on the temperature. 
At Tf r w 0.8 MeV, the n/p ratio freezes out and thus 

(— J « exp(-Q/T fr ) « 0.2. (2.20) 

More exact calculations || give (n n /n p ) {r « 1/6. 

After freeze-out, the n/p ratio continues to decrease slowly due to neutron 

decay, 

n^p + e"+z/ e . (2.21) 

Therefore, 



, / BBN 




ln(2) tBBN 



~, (2-22) 



where r n = 886.7 ± 1.9 s is the neutron half-life. If we insert this value into 
d2l7D , we get Y w 0.25. 



2.4 Non-Standard Neutrinos 

We now want to derive the influence of neutrinos on the helium abundance. 
According to (|2.19|) , a non-zero electron neutrino asymmetry changes the 
n/p ratio by a factor 

exp(-/v/T fr ) w exp(-1.5^J w (1 - 1.5^J (2.23) 

provided that A Uc < 1. As an example, we take A Vc = ±0.01, a realistic 



value according to Ul9[ . Then (n n /n p ) is altered by a factor (1 =f 0.015). Of 
course, this is a very rough estimate, but more detailed works JTB], |37J find 
the same order of magnitude. 

A higher effective neutrino number also alters (n n /n p ), since it raises g*, 
which results in a higher expansion rate H and thus in a higher freeze-out 
temperature. Taking = 4 instead of 3 as an example, Tf r is increased by 
a factor of 



' g ,(N? = 4) \ 1/6 = / 10.75 + 2 • | \ 1/6 

M N f = *)J \ 10 - 75 / 



1.025, (2.24) 
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where the power 1/6 comes from Y/H oc T 3 / y/gl. This changes the n/p ratio 
by a factor of 1.04. 

We insert these changes into ( |2.17[ ) and expand to get 



Y w 0.25 + 0.01(iV^ ff - 3) - 0.33A, o (2.25) 

if A Ve < 1 and \Nf - 3| & 1. We see that A^ ff = 4 or |A/ e | = 10" 2 change 
the helium abundance by several percent, which is within the present exper- 
imental precision. More detailed numerical calculations || give 

Y w 0.225 + 0.025 log^/lO" 10 ) + 0.012(^ ff - 3), (2.26) 

where we have also included the influence of the baryon asymmetry 7], since 
the main uncertainty in Standard BBN comes from this parameter. 

We should mention that in the literature, the variable N* s has often been 
used not only to account for the effective number of neutrinos, but also for 
the influence of A Vc . This was done by adding a new term SN^^A^) to Nff . 



2.5 Observational Results 

The best measurements of the primordial 4 He abundance, Y, come from 
observing extra-galactic regions of ionized H. In these systems, the abundance 
of heavier elements, which are not created in BBN, is very low, so we can 
assume that the abundances in these regions are close to their primordial 
values. The present estimate is 

Y = 0.238 ± 0.002 ± 0.005, (2.27) 

where the two errors are the statistical and systematic errors, respectively. 
The 2cx range is then estimated @ to be 0.228-0.248. 

Direct present-time measurements of the cosmic baryon abundance are 
very uncertain due to the dark matter problem. Much better estimates 
arise from measurements of the primordial Deuterium abundance which de- 
pends sensitively on the baryon asymmetry rj. The best measurements come 
from the absorption of quasar light by high-redshift, low-metallicity hydrogen 
clouds. Two main results have been published [^, ffij, 

(D/H)i ow = (3.4 ± 0.3) x 10~ 5 and (D/H) hi gh = (2 ± 0.5) x 10" 4 , (2.28) 
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which are mutually inconsistent. These two measurements give ranges for 
the baryon asymmetry of 



Low 



= 4.2-6.3 x 1(T 10 and r] high = 1.2-2.8 x 1(T 10 , (2.29) 



at a nominal 2a level. 

From the measurements of D/H and Y, one obtains bounds on the effec- 
tive neutrino number. Olive et al. [fUl derived 



(Nf - 3) low < 0.3 and (N? - 3) high < 1.8, (2.30) 

provided that the v e chemical potential can be neglected. 

The low-D result today appears to be strongly favored. If it should be 
confirmed, sterile neutrinos would be forbidden to come into equilibrium for 
negligible fi Uc . However, a large positive A Ve can compensate the effect of an 
increased iV7; ff , and can thus circumvent this constraint. 



2.6 Summary 

We have shown that a deviation of the effective number of neutrinos of order 
1, as well asymmetry exceeding about 10 2 , will have a measurable 

effect on the outcome of BBN. Sterile neutrinos have the potential to change 
both Nf 1 and A Uc . If we knew how these two variables depend on the mixing 
parameters of the sterile neutrinos, we could use the measured primordial 
element abundances to derive constraints on the mixing parameters. On 
the other hand, if the existence of sterile neutrinos was proven by future 



experiments such as MiniBooNE |]4"2" |, a detailed understanding of their im- 
pact on primordial nucleosynthesis would be necessary to constrain the free 
parameters of BBN, in particular the baryon asymmetry rj. 
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Chapter 3 

Neutrino Oscillations 



Neutrino oscillations play a crucial role in our considerations of the Early 
Universe involving sterile neutrinos. In this chapter, we derive the density 
matrix formalism for neutrino oscillations in media between any two neu- 
trino flavors, active or sterile. This includes the direct influence of matter on 
the vacuum neutrino oscillations as well as scattering processes which tend 
to destroy the coherence of the oscillations. Our analysis will only be ap- 
plicable for temperatures between a few MeV and 100 MeV and for neutrino 
asymmetries A v <C 1. 

3.1 Equation of Motion 

Neutrino oscillations occur because the basis of the neutrino weak eigen- 
states v a , a = e, /i, r, s, . . ., is different from the basis of the neutrino mass 
eigenstates i = 1,2,3,.... In other words, a v a that is produced in a 
weak-interaction process does not propagate like a free particle, but as a su- 
perposition of neutrinos vi with different masses raj, respectively. Thus, when 
measured, the propagated neutrino contains contributions of weak eigenstates 
other than the original v a . From the start we include the possibility that neu- 
trinos exist beyond the active states u e , and u T . These additional flavors 
would have to be sterile with regard to the weak interaction. 

The two bases are connected by a unitary transformation \l/ w = U^m, 
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where 



w 



and * 



M 



"2 



(3.1) 



V 



/ 



V 



/ 



are the field vectors represented in the weak and mass basis, respectively, 
and U is the unitary transformation matrix. 

We are mainly interested in the evolution of the weak eigenstates, as 
these are the particles we can produce and measure by weak-interaction pro- 
cesses. Therefore, we write the equation of motion in the basis of the weak 
eigenstates. The Klein-Gordon equation is 



( 5 2_ V 2 + M^)^ = 0, 



(3.2) 



where Mw = UMyJJ^ and Mm = diag(rai, m-z, ■ ■ •) is the mass matrix in 
the weak and mass basis, respectively. Of course, = UMuU^UMuU^ = 
UMyJJ* and M\y are n °t diagonal in the case of mixing. 

Since the Early Universe is homogeneous, we are only interested in the 
time evolution of the neutrino fields. It is therefore convenient to expand 
them in plane waves \l/w — ^ p (t)e ipx . Note that usually neutrino oscil- 
lations are considered in environments which are spatially inhomogeneous 
but stationary, e.g. in experiments involving solar or atmospheric neutri- 
nos. In such cases, \&w would be expanded in components of fixed energy, 
\l/w — \pE(x)e~* , instead of components of fixed momentum. 

Since the neutrinos are highly relativistic, m v <C E ~ p, we can linearize 
the Klein-Gordon equation. Then we get the usual Schrodinger-type equation 



= V + 



2p ' 



(3.3) 



We see that the off-diagonal elements of M(y couple the fields of the different 
neutrino flavors, leading to oscillations. 
An equivalent equation is given by 



id t p = [p,tt], 



(3.4) 



where p aj 3 = N^a^p is the flavor density matrix, N is a normalization factor, 
and we have dropped the index p. The advantage of the density matrix 
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formalism is that we later can include effects that destroy the coherence of 
the neutrino oscillations. 

For later convenience, we normalize the density matrix such that for a 
CP symmetric neutrino species u a in equilibrium, p aa = 1 for all momentum 
modes, i.e. 

where fo(p,T) = f(p,T,fj, = 0) = [1 + exp(p/T)] _1 . Thus, the normalization 
factor N = l//„(p,T). 

We have neglected the expansion of the Universe in the derivation of the 
equation of motion. It could be included by adding a term fl H = iHpd p 
to the Hamiltonian. The term disappears if we expand the field vector in 
comoving plane waves, = ^ij(i)e 8Tqx , where q = p/T, instead of plane 
waves with fixed momentum. Therefore, later in Chapter 4 we will not 
describe the evolution of flavor density matrices of fixed momentum p, but 
of fixed comoving momentum q. 



3.2 Medium Effects 

In the Early Universe, we are confronted with a thermal medium which 
interacts with the neutrinos. Therefore, we need to include the effects of 
these interactions on the neutrino oscillations. In this section, we restrict 



ourselves to the discussion of the refractive effects [f45|j . 

The medium contributions to the neutrino oscillations enter the Schrodin- 
ger equation ( |3.3| ) through the weak-potential term V = di&g(V Ua , V Uf3 , . . .) 
in 

M 2 

n = p+-™+V. (3.6) 
2p 

For each neutrino weak eigenstate the contributions can be split into two 
terms, 

v Va = ±v£-v£. (3.7) 

Here, the plus sign is valid for neutrinos, while the minus sign applies to 
antineutrinos. 

The first term accounts for the fermion asymmetries. For neutrinos of 
species a it is given by 

V* = V2G F n,A Uci , (3.8) 
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where Gp = 1.166 x 10 5 GeV 2 is Fermi's constant and A Va a weighted sum 
over all fermion asymmetries Aj. Generally, 

A Va = A Va + A Vc + A U)i + A„ T 

+ A a - -(1 - 4 sin 2 0w) (A e + A M + v4 T ) 

+ l(l-4sin 2 M^ P -^n, (3.9) 

where #w is the Weinberg angle. For our assumption of a charge-neutral 
Universe we have A e + A^ + A T = A p and therefore 

A Ua = A Ua + A Ve + A Vfi + A Vt - ^A n + A a . (3.10) 

Since the muons and tauons are non-relativistic, their asymmetries are neg- 
ligible, so that we will use A e = A p , A^ = 0, and A T = 0. In our analysis, 
all of the asymmetries but for one type of neutrino will remain constant. 
Therefore, we write 

A Va = 2A Va +A c , (3.11) 

where all asymmetries other than A Va have been absorbed in a constant A c . 

The second term in ( |3.7|) represents the low-energy tail of the W and Z° 
resonances. It has the remarkable feature that it is independent of the CP 
asymmetry of the background medium and has the same effect on neutrinos 
and antineutrinos. Since it is of order Gp, V T can only compete with V A in 
the Early Universe, where A<C 1. For neutrino temperatures far below the 
W ± and Z°-masses it can be written in the form f45 



Here, the coefficients are 



Vl = z„ a G 2 F pTn T (3.12) 



77T 3 

Zu, = sin 2 # w (cos 2 # w + 2) w 350 

45C3« v 7 

7vr 3 

z v „, Vt = —— — sin 2 6\y cos 2 ~ 97, (3.13) 
45C 3 a 

where a = 1/137 is the fine-structure constant. 

Since we will be discussing sterile neutrinos, we stress their special role in 
a medium. They do not interact, so V Us = 0. Likewise, a possible asymmetry 
A,, does not contribute to V A for the active flavor. 
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3.3 Density Matrix Formalism for 
Two-Flavor Oscillations 

We will now treat the simplest case of neutrino oscillations between two types 
of neutrinos, so — {^v a i ^v )- The unitary transformation can be written 
in the simple form 

TT ( cos^ — sin# \ 

U = • a a ' 3 - 14 
\ suit/ COS# / 

where # is the vacuum mixing angle. Note that the absolute phase of the off- 
diagonal terms is arbitrary; our convention differs from some of the literature. 
In the weak basis we get 

il=-bI--( cA+ . SV , A SA „A (3.15) 
2 2\sA -(cA + 5V) r y ' 



where b = 2p + mi 2™ 2 — (K>„ + K/g), I is the 2x2 unit matrix, and 5V = 
V Va — V U0 . Furthermore, we have defined 

Styi 2 

s = sin2#, c = cos2^, A = (3.16) 

2p 

with the mass squared difference 5m 2 = m\ — m\. 

To find a convenient form for the equation of motion, we represent p and 
Q by the unit matrix I and the Pauli matrices cij, i = 1, 2, 3 [PJ]. Then 

P=\P Q I+\?(T, n = ±bI-^B<T, (3.17) 
where P is a so-called flavor polarization vector, and 

B = | I . (3.18) 




Inserting p and f2 into fl3.4p and using the commutation relation for the Pauli 
matrices, [|<7j, |<x,] = ie^k^Uk, we derive 

<9 t P = BxP. (3.19) 

This equation is equivalent to the precession of a spin vector in a magnetic 
field. The physical observables are now contained in P z = p aa — ppp and 
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Pq = Paa + P/3/3- fo(p,T)Po is the total occupation number of both v a and 
vp of momentum mode p, while P z represents the difference between the two 
neutrino types; if e.g. P z = 1 for all modes, then n Vfj = 0. However, if P z = 
for all modes, n Va = n Vf3 . The other two components P XyV represent the phase 
of the oscillations. 



3.4 Scattering Processes 

We must still take into account that the weak neutrino eigenstates scatter 
on the background plasma with the rate |46| 



r, Q = y Va G 2 FP T\ (3.20) 

where y Uc ~ 1.13, y Vll ,v T ~ 0.79 and y Vs = 0. To be precise, y Va also depends 
on the momentum due to the Pauli blocking factors. We will use these 



momentum averaged values [24 



Scattering keeps the active neutrinos in thermal and chemical equilib- 
rium at temperatures above a few MeV. Therefore, the integrand in ( p.5|) , 
fo(p, T)p aa (p), should be equal to f(p, T, /i). However, p a a(p) changes due to 
the neutrino oscillations. The scattering processes compensate this change 
by re- or depopulating the number densities through p aa (p) = R Ua (p), where 
the rate is approximately given by p6f 



\(p)«r t 



f{p,T,p) , , 

~ Paa{P) 



(3.21) 



fo(p,T) 

The approximation holds as long as the deviation from equilibrium is small. 
Of course, R Us = T Us = 0. 

A much more interesting effect appears if the scattering rates are dif- 
ferent for the two neutrino types considered. Then scattering distinguishes 
between the two v types and thus the coherence of the neutrino oscilla- 
tions is destroyed. Therefore, we need to include a term —DP±(p), where 
Pi = (P x , P y , 0) and D is called the damping rate. In the case of an active 
neutrino oscillating with a sterile neutrino ((3 = s) we have D = \ T Ua . The 
lengthy derivation of this term can be read in |[47|| . 

In summary, for each neutrino momentum mode p we have derived a 
system of differential equations 

dtP = B x P — DP± + (R Ua — R U/3 )z, (3.22) 
d t P = (R Ua +R^), (3.23) 
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d t P = B x P - DP ± + {R Va - R Vf> )z, (3.24) 
d t P = (R Ua + R Up ), (3.25) 

where z = (0, 0, 1). In the next chapter we will study the highly non-trivial 
behavior of this equation of motion. 
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Chapter 4 

Oscillating Lepton Asymmetry 



We analyze the system of flavor oscillations between active and sterile neutri- 
nos before BBN in a simplified model. First we discuss our approximations 
and present a typical numerical solution of the simplified model. Next, we 
introduce a more convenient coordinate system for the analytical treatment. 
With its help, we describe the evolution of our system from very early times 
up to the resonance. At resonance, we show that we can describe the evolu- 
tion of the neutrino asymmetry with a simple differential equation, and that 
the solution will indeed oscillate for some of the parameter space. 



4.1 Simplified Model 

We describe the oscillations between two flavors of neutrinos, one active 
and one sterile, in an expanding medium. This means that we neglect all 
other neutrino mixings, which is a good approximation if all other effective 
mixing angles are small. Thus, we can describe our oscillations by the system 
of differential equations ( |3.22j )-( |3.25| ) derived in the previous chapter. We 



restrict ourselves to a simplified model which we present in this section. For 
definiteness, we will analyze the case where the active neutrino is v T \ the 
analysis is analogous for u e and v^. 

Our most important simplification is that we neglect the momentum dis- 
tribution. This is surely not a good approximation, since then all neutrinos 
encounter the oscillation resonance at the same time. In reality only a small 
fraction of the neutrinos will be simultaneously at resonance, especially when 
the resonance width is small, i.e. the vacuum mixing angle is small. However, 



23 



the complete system is very complicated, as can be seen from the controver- 
sial literature on its solution, and we therefore have decided to analyze this 
simplified model. Of course, it is desirable to include the effects of the mo- 
mentum distribution in future investigations. 

So from now on, all neutrinos are taken to have the same momentum. We 
choose this momentum to be the average of a fermion species with vanishing 
chemical potential, 

<p) = ^T^3.15T, (4.1) 

where ( is the Riemann zeta function with £ 3 « 1.202 and £4 = 7r 4 /90. As 
a consequence, we only have two density matrices, p aa for the neutrinos and 
p aa for the antineutrinos. We normalize these in analogy to our previous 



normalization in Section |3.1| , i.e. such that the number densities are n Ur = 
nl q p TT , etc. where n eq = |n 7 is the equilibrium neutrino number density for 
vanishing chemical potential. 

We will also neglect the repopulation terms R Ut and R Ut defined in ( p.21|) . 
They will be small if the v T and v T are not depopulated significantly by 
neutrino oscillations, and if A Vt C 1 so that the chemical potential has no 
significant influence on the equilibrium number density. This simplification 
is valid for some part of the parameter space, which we will determine later. 
As a result of this approximation, Pq and Pq are constant in time. 

So now, the system of differential equations has been simplified to 

dtP = BxP-DPi, 

d t P = B x P — DP±. (4.2) 
The tau neutrino asymmetry 



eq 



n 

A Vr = -^(p rr - p TT ) = I x UP Z -P z + P - p Q ) (4.3) 
n 7 8 2 

is an important variable since it depends on P, enters into B, and thus causes 
the system to be nonlinear. It will be more convenient to use the equivalent 
variable 

6P z =±A Ur = ±(P z -P z )+P G , (4.4) 



where A Vt was given in ( |3.10| ) and P c = §(A/ e + ^ M _ + h(Po - Pq) is 



2^ a > 1 2' 

constant. 
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In terms of this variable, the coefficients in the differential equations are 
given by 



B(T,6P Z 
B(T,5P Z 





-b c + b T - b A SP z 

-b s 


V -b c + b T + b A SP z 



(4.5) 
(4.6) 



where we use 



s 5m 2 
c 5m 2 



s 5m 2 

6.3 
c 5m 2 



2(p) 6.3 
b A = V u A J5P z = k A T 3 , 



b T = Vl = k T T\ 
Recall that s = sin 28, c = COS26 1 , and 5m 2 



m 



(4.7) 

1 — m 2 . The constants are 



k A 
k T 



3(3 

4tt 2 
14vr 

45a 



1.51 x 10" 24 eV 



sin 2 6 W cos 2 6 W ^-G 2 F « 1.02 x HT^eV -4 , 



(4.8) 



where a = 1/137 is the fine-structure constant and 6\y the Weinberg angle. 
The coefficient D is proportional to b?, so we will often use the relation 



D = k D b T , 



where 



1/60, 



(4.9) 



(4.10) 



and we have used y Ur ~ 0.79 and z Vt « 97 from Section |3~2"[ For u^, the 
constants are the same. For u e , we have k^ ~ 3.66 x 10 _44 eV -4 and ko « 
1/151. 

We will be considering neutrino oscillations with small vacuum mixing 
angles, i.e. s < c « 1. As a consequence, the first component of B will be 
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much smaller than the third component for most of the time. The system be- 
comes interesting when the third component disappears. Then the neutrino 
oscillations are at resonance, i.e. the neutrinos mix maximally. If we assume 
the initial asymmetry to be negligible, the resonance condition is given by 
b c = br- This condition has a solution only if 5m 2 < 0. Then resonance 
occurs when 

T res = (^) V6 « 15.8MeV|5rr4r, (4.11) 

where 5m 2 Y = 8m 2 / eV 2 . The resonance will be the crucial feature of the 
system. We will see that shortly after the resonance, the neutrino oscillations 
will create a large asymmetry, an effect which is driven by the non-linear term 
bA5P z in the system of differential equations. We will therefore only consider 
Sm 2 < 0, i.e. the sterile neutrino is lighter than the tau neutrino. 



4.2 Numerical Solution 



We have solved the system (4.2) of differential equations numerically and 
find results similar to those in and ||22|| . We have plotted the evolution 
of the effective asymmetry 5P Z for (5m 2 , s) = (— leV 2 ,10 -4 ) in Fig. 4~T 



Its behavior is representative for a large region of parameter space. The 
evolution of the system falls into five distinct phases. 

1. At high temperatures, the asymmetry 5P Z is constant, since damping 
totally destroys the coherence of the oscillations. Thus, no flavor tran- 
sition occurs at all. 

2. As the temperature decreases, the mixing angle increases, while the 
damping rate decreases. As a consequence, there will be a small amount 
of flavor transition induced by the damping. Due to the asymmetric 
term b^SPz, this transition rate will be different for neutrinos and an- 
tineutrinos, washing out the asymmetry 5P Z . However, there remains 
a small relic asymmetry. 

3. Until the system reaches the resonance, the relic asymmetry changes 
slowly. 

4. When the system has passed the resonance, the difference in the tran- 
sition rates for neutrinos and antineutrinos has the opposite effect than 
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logio(l^l) 




Figure 4.1: Evolution of 5P Z , steps 1. to 4. 



<5P,[10- 7 ] 




15.78 15.76 15.74 15.72 



15.7 



T[MeV] 

Figure 4.2: Oscillatory behavior of SP Z during step 5. 
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before. In other words, a small asymmetry is no longer washed out, 
but amplified, which leads to an exponential increase of 5P Z . 

5. When SP Z = 0(1O~ 6 ), it starts oscillating and thereby even changes its 
sign. For our numerical example, we have plotted the evolution of 5P Z 
during this step in Fig. 



6. Eventually, the neutrino asymmetry stops oscillating. After that, it 
slowly increases. 

In the past, the underlying mechanism for the oscillatory behavior in 
step 5 was not fully understood. Merely the authors of |21| discussed this 



phenomenon qualitatively. Our main achievement will be to analyze the 
system analytically and to prove that it indeed oscillates. Thus, we invalidate 
the frequently used argument that the oscillations in this phase are an artifact 
caused by a numerical instability. 



4.3 Spherical Coordinates 

We will now change to a more convenient coordinate system. The two phase- 
related coordinates P x and P y are not very practical, so we prefer to use 
spherical coordinates. As a measure for the length P = |P| of the polarization 
vector we use P z because it enters directly into ( |4.4| ). This choice becomes 
problematic when P z = 0, P^O, but as mentioned above, P z will not change 
significantly. The new coordinate system consists of P z and the angles $ and 
ip. These new variables are related to the old ones as 

P x = P z tan sin <p, 
Py = -P 2 tani?cosv3, 

Pz = Pz- (4-12) 

P z represents the relation between the number densities of the active and 
sterile neutrinos, °d is a measure of how much the polarization vector deviates 
from the z-axis, and ip is the rotation angle around the z-axis. Then the 
differential equations (|4.2j ) change to 

P z = —b s P z tan $ cos <p, 
*d = — D sin^cosi? + b s cosip, 

(p = b c — bx + bA^Pz — b s cot $ sin <p, (4-13) 
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and 



P z = —b s P z tan •& cos <p, 
•d = —D sin •§ cost? + b s cos (p, 
<p = b c — bx — bj^SPz — b s cot d sin (p, 



(4.14) 



we 



where over-barred quantities, as usual, refer to antineutrinos. With 
have a closed set of differential equations. 

The only difference between the equations for neutrinos and antineutrinos 
is the sign of the term bA$P z - Therefore, if this term is negligible, the variables 
for neutrinos and antineutrinos develop equally. To be more sensitive to the 
small differences induced by the asymmetry, we perform another substitution 
and use the variables 



Pf 



instead. Thus P+ 3> P~ 



$ + ^> $~ and ip + ^> ip~ when the fermion asym- 



(4.15) 



metry is small. Note that Q4.4 ) becomes 

5P z (t) = P z -(t) + P c 



(4.16) 



so that ^8P Z = P~ . We will make use of this simple relation by considering 
5P Z instead of P~ when it is convenient. 

In terms of the new variables, our system of differential equations is 



p: 



-kp: 



-b s P z ~ 



-b P + 

u s 1 z 



-b s P, 
1 



sin 2$ + 


cos <p + cos ip 


— sin 2$ 


sin ip + 


sin (p~ 




cos 2t?+ 


+ cos 2tf~ 






-sin2$ H 


" simp + sin<y?" 


" + sin 2d~ 


coscp + 


cos ip~ 




cos 2d+ 


+ cos 20- 






-sin 2^ 


' sin (p + sin (p~ 


" + sin 2d~ 


cos<^ + 


cos (p~ 




cos 2$+ 


+ cos 2^" 






sin 2$+ 


cos (p + cos ip 


" — sin 2d~ 


sin ip + 


sin (p~ 



cos 2$+ + cos 2$" 
D sin 2d + cos 2-$~ + b s cos ip + cos 



2 

--Pcos 2$ + sm2iT 
2 



6 S sin ip + sin 93 



(4.17) 

(4.18) 
(4.19) 
(4.20) 
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cos 2i9~ — cos 2-$+ 

(p~ = b A {P; + P c ) 

— sin 2d + cos 09 + sin tp~ + sin 2d~ sin 09 + cos 09" 



+ b s 



cos 2i9~ — cos 2§ + 

(4.22) 



As a last approximation, we assume d ± <C 1, which effectively means that 
the polarization vectors will stay close to the z-axis. We expand sin 2tf± and 
cos 2$ ± and only take into account the leading order terms. Then we get 



= —b s P^ ( $ + cos 09 + cos (p~ — d~ sin 09 + sin <p~ 

—b s P~ (— d + sin 09 + sin 09" + i?~ cos 09 + cos o9~) , (4.23) 

P~ = —b s P^ (— $ + sin 09 + sin o?~ + *&~ cos 09 + cos o9~) 

—b s P~( $ + cos 09 + coso9~ — "&~ sin 09 + sin o9~) , (4.24) 

= + & s cos o9 + cos y9~, (4-25) 

0- = -Di?" -6 s sinv9 + sinv9", (4.26) 



+ ^ ( - sin 09 + cos ip + — cos y9 + sin 09 j, (4.27) 



<j>~ = b A {p- + P c ) 

b s ( , _ iT , _\ , 

+ —— cos 09 sin 09 + — sin 09 cos 09 . (4.28) 
y ' 17+ y 

From (|Q5D we see that < -D&+ + b s . Thus i?+ > b a /D would auto- 
matically mean that $ + < 0. Therefore, we can safely say that i9 + < b s /D 
at all times. We will see that $ + becomes maximal at resonance, so ~ 
(b s /D) res = s/(ck D ), where we have used b s /s = b c /c, (6 c ) rcs = (b T ) res and 
D = kobr- Thus our approximation is only valid if s -C ko = 1/60, c ~ 1, 
i.e. for small mixing angles. 
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4.4 Initial Conditions 



We need to derive the initial conditions before we can consider the evolution 
of the system of differential equations. We begin with temperatures far above 
the resonance. For T — ► oo, the coefficients D, br, R Ur oc T 5 — > oo, 6^ oc 
T 3 — > oo, while b s , b c oc T _1 — > 0. If we compare these coefficients with H oc 
T 2 , we conclude that b s and b c can be neglected at very high temperatures. 
Thus, from ( 4.23|) - (|4.26|) we see that regardless of the initial conditions, $ + 



and $~ are both exponentially damped to zero, while Pf oc b s do not change. 
However, the scattering processes equilibrate v T and z/ T , implying p TT and 
p TT — > 1. Since the sterile neutrinos do not interact, p ss and p ss remain 
constant. They are only diluted whenever massive particles become non- 
relativistic and annihilate into the still relativistic particles, which heats up 
the plasma. For example, when the temperature decreases from T ~ TeV 
(just before the electro-weak symmetry breaking) to T ~ MeV (shortly after 
the quark-hadron phase transition), p ss and p ss are diluted by a factor of 
[#*(T ~ 300MeV)/#*(T ~ TeV)] 4 / 3 ~ 0.07. So if the u s were in equilibrium 
with the thermal plasma at very early times by some unknown mechanism, 
they would have been strongly diluted by the time they become interesting 
for us. We will simply assume that their initial density is zero. Then the 
initial conditions are Pq = Pq = 1 and P z — P z — 1. Including the small 
initial neutrino asymmetries, A Vt and A Us , we find 

p+ = l 
4 

P z = -jv^v — A Va )i 
1-4 

- (P - Pq) = -(A Ut +A Us )= const , 

P c = ^{A c + A Ut + A Us ) = const, (4.29) 
where in P+ we have neglected A Vr and A Vs <^ 1. 



4.5 Quasi-Static Solutions 

We will now follow the time evolution of the system. We first stress the 
important role that damping plays in our system. We consider a variable x, 
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which follows the differential equation 



l x (t) = -d(t)x(t)+f(t), (4.30) 

where d(t) is some damping coefficient and f(t) is a function. If d and / are 
constant, x will relax to a static value x st = 4. In our case, d and / will slowly 
vary in time, and so will the static value. But if its rate of change r x = x st /x st 
is smaller than the damping coefficient d, the damping will force the variable 
to follow its static value. We will make use of this approximation; since it 
is not static in the strict sense of the meaning, we will call these solutions 
quasi-static. 

We have calculated the rates of change of the quasi-static solutions for 
our variables in Appendix For the discussion in this section it is sufficient 
to know that r = 0(H). 

We can now determine the variables for large temperatures. We can easily 
assume that <C Dfl^, b s , (b c — br),^/^, so before Pf change signifi- 
cantly at all, the other four variables will have relaxed to their quasi-static 
values. If we assume that P~ = 0, the differential equations for neutrinos 
and antineutrinos ( |4.2j ) decouple. For small the angular equations in ( |4.13| ) 
become 

= -D& + b a cosy (4.31) 

y = (b c -b T + b A 6P z )-^sm<p. (4.32) 

v 

Since d > by definition, the first term in the first equation and the second 
term in the second equation are both damping terms. The two variables 
relax to « ip if the quasi-static conditions are fulfilled, i.e. D > H 
and bg/d > H, respectively. The first condition is fulfilled for T > fewMeV. 
When t? has relaxed, we get from (|4.31| ) that ~ b s cos f/D < b s /D, which 
fulfills the second condition, b s /$ > D > H. Setting $ = (p = 0, we get the 
quasi-static solutions 

tf qs = : bs (4.33) 

Id* + (b c - b T + bjpr 



f b c -b T + b A 5P z . 
(p qs = arctan I — ) . (4.34) 
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The same is valid for the antineutrinos, but with 5P Z — > —5P Z . Since 
b A 8P z "C b c — bx, we can already say that the variables $~ and <p~ will 
be much smaller than i9 + and (p + , respectively. Thus, 

#t = |(^ qs +^ qs )~ b ! = , (4.35) 



qs 



-D 2 + (6 C - b T ) 



qs 



|(^qs + w) ~ arctan 



~ b T 



(4.36) 



The same quasi-static solutions are also derived from ( 4.25 ) and (|4.27|) if we 
neglect the last term in latter. 

Now it is also easy to calculate $~ and <p~ for a given SP Z by setting 
$~ = in ( |4.26| ) and ip~ = in ( |4.28| ), inserting i?+ and and setting 
cos(/?~ = 1. Then 



D6, 



<^qs 



arcsm 



D 2 + ( &c _ &T )2 

b s (b c -b T )b A , n 
did. 



(4.37) 

[ J D 2 + (& C -M 2 ] 3/2 ^ 2 ' (4 ' 38) 

Again, we need to show that the quasi-static conditions are fulfilled, which 
we again find to be D > H. 

If we insert all these variables into P~ , and neglect the fourth term in 
( |4.24j ), we obtain 



d_ 

dt 



5P z = P- = -K{5P z + eP-), 



where the damping coefficient is 

2P+b 2 s (b c -b T )Db A 



K 



and 



[D 2 + (p c - b T f} 2 
D 2 + (6 C -6 T ) 2 



oc T 



-9 



" 2P+b A (b c - b T ) 
where e <C 1 for T < GeV and T 7^ T res . Analogously, we get 



(4.39) 

(4.40) 
(4.41) 
(4.42) 



where we have only taken into account the first term in Q4.23 ), since the 
other terms are all 0(P~) smaller. This confirms that P~ and P z will not 
relax to their quasi-static values at high temperatures. 
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4.6 Evolution Towards the Resonance 

When k^H, 5P Z starts changing. The damping coefficients for ip~ and i? _ 
(which are ^> D) are much larger than H, so even though 5P Z changes, (|4.37|) 
and (|4.38|) remain valid. 

From ( 4.39|) we see that 5P Z is damped towards 0, i.e. P~ — > — P c . When 
\SP Z \ <C \Pz\i the second term in ( |4.39| ) becomes important. We set P~ = 
to get 

(p-) qs = -P c (5P z ) qs = P c (4.43) 

We have to prove that SP Z takes on its quasi-static value before the system 
passes the resonance. The solution of ( |4.39| ) is given by 

5P z (t) = 5P Z (0) exp ^- J re(t)cftj . (4.44) 

If we substitute x = T(t)/T res , we can write the integral as 



where 



,211/6 



6.3V8 5.5 k T 



Assuming that initially 5P Z and P~ are of the same order of magnitude, we 
know that 5P Z will decrease by a factor of order e ~ &t/^a ~ 0(1O~ 6 ), which 
corresponds to P_(t res ) > 61n(10). Thus, we obtain the condition 



S' 



'ItfrnV^^lO-^eV 1 / 3 (4.47) 



using P+ « 1. 

For P+ we can perform a similar calculation. However, here we demand 
that P 2 + does not change significantly to prevent that the repopulation terms 
become important. Therefore we demand that F + (t res ) <C 1, where 

P/(t) = P+(0)exp[-P+(t)]. (4.48) 
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Thus, 



where 



'l /2 «2|A m 2|V2 



/V(/) i^t)dt = !^^ s -^-un. tun) 







;r 2 



I+fe.) = / k>12 - (g6 _ 1)2 rfx « 15.2. (4.50) 



We then get the condition 

s 4 \5m 2 eY \ «2x 1(T 9 , (4.51) 

where Sm 2 ^ = 8m 2 / eV 2 . F + (t ves ) < 1 is also the condition that the v B do not 
come into equilibrium. Our result here is in good agreement with previous 
constraints [IS . 



The quasi-static approximation is of course not exact, actually the vari- 
ables will always be a bit behind their quasi-static values. This delay becomes 
important when the quasi-static value of a variable passes zero, or when the 
rate of change of the quasi-static value becomes of order of the damping rate 
of the variable. At resonance, i.e. when (b c — fey) = 0, we have d~ s = 0, 
V^qs = 0' (Pz)qs — 0. Furthermore, all the quasi-static values change rela- 
tively fast near the resonance. 

Here is a qualitative discussion of the behavior of the variables close to 
resonance. The term that is responsible for the resonance, Afe = b c — fey, 
only enters into the differential equation of ip + . When Afe changes sign, y?+ 
does too. Therefore, after a short delay, ip + will also change its sign. The 
other variables do not depend directly on Afe; they only feel its change via y? + . 
Thus, they are not sensitive to the time delay between (p + and tp~L. Therefore, 
setting ip + = (fi+ s will merely shift the events by a negligible amount of time. 

Close to resonance, we have sinip + w and cos(/? + ~ 1. Therefore, $ + 
becomes maximal close to the resonance, i.e. $ + ~ b s /D, while $~ becomes 
minimal so that we can neglect it in the equations. tp~ will stay close to 
its quasi-static value since 5P Z does not change significantly: SP Z freezes out 
when the rate of change in (SP z ) qs becomes larger than the damping rate k, 
i.e. rsp z > k, where r$p z is given in ( A.7Q . Since Afe = fe c — 6 T <ti D and 
e <^ 1 at freeze-out, we get 

nT b c 2P+b 2 s b A Ab , Arn . 

r 5Pz « -6#^, k w 2 ^ 3 • (4.52) 
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Using T = T res , we find that Ab h = -8.1 x lO-^-^m 2 ^- 1 / 12 ^ 1 / 12 ^, and 
thus e b = 7.2 x l(r 3 s|<5m 2 v | 5 / 12 c- 7/12 < 1, so that (6P z ) b « e h P c is indeed 
very small. 

In summary, at resonance we can still approximate the variables $ + , <£> + 
and ip~ by their static values, while we can neglect 5P 2 is given by its 
freeze-out value. 



4.7 Evolution at Resonance 

We now come to the most interesting feature of the system: sin (p + changes 
its sign and becomes positive. Therefore, the back- reaction of (p~ on 5P Z , 
which is given by the coefficients of the first terms in ( 4.24 ) and (|4.28 ), 



bA • b s P^'d + sin also changes its sign to positive. So while 5P Z until 
now was washed out through this back-reaction, the same back-reaction now 
amplifies SP Z . We will now derive the consequences of this effect. To this 
end, we will assume that the variables ^ and ip^ will be given approximately 
by their static values. As time passes, these approximations will break down 
and we will need to find other approximations. 

4.7.1 Static Approximation after Resonance 

The change of 5P Z becomes interesting again when r$p z < k after the reso- 
nance, i.e. when Abo — —Ab{ r . We will denote the time when this happens 



with t - So until then, 5P z (to) ~ (5P z )f r = S{ r P c . Now (|4.39|) has changed to 



j t 5P z = +k(5P z - eP~) « +k(5P z + eP c ). (4.53) 
If we neglect the second term, we get the solution 

5P z (t) = 5P z (t ) exp |y" Kdtj . (4.54) 

To solve the integral, we substitute t with £ = 1—x <^ 1, where x = T(t)/T Tes . 
Then 

Ab = b c -b T = fe^aT 1 - b^ s x 5 = b^ s x-\l - x 6 ) w 66?r s £, (4.55) 
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and 

d£ = --^dT = ^^ dt « ^^dt, (4.56) 
-/res m p i m p i 

where we have expanded £ and have only taken into account the leading 
order. The integral becomes 

\dt = J^-/ gg%ggi de 

5.5T2 7 ^ 5.5T2 7 £3 ^ 



, ■ res r res ,. 

*0 fo 5o 



2P+b 2 s b A 6b T m p i 



D 3 5.5T 2 

5.5 x 10 17 S 2 |^ev| 1/6 c- 11/6 • |(^ 2 - ( 4 - 57 ) 



T-T rcs ?o 



where we have assumed that Ab -C D and have again expanded £ to leading 
order. Furthermore, £ = A6 /66^ s = 1.35 x 10 _9 s~ 1 |5mg V | _1 / 12 c 11 / 12 . 

The solution ( J4.54|) is only valid as long as <p~ and i9~ can follow their 
static approximations, which are given by ( 4.37|) and (|4.38|) . We see that 



as long as Ab = b c — br <C D, the changes in Lp qs and $ will mainly be 
due to 5P Z , i.e. (p~J(p~ s w ^qs/^qs ~ &P t /6P g . Then the validity of (p§ 
breaks down when D = -^5P Z /5P Z = k, which happens when £ = £i = 
1.1 x 10 _14 s _2 |5m 2 v | 1//3 c 7//3 . We see that the above solution already breaks 
down very early, for a large parameter range we even have £i < £o- To give 
an estimate of how big the integral ((4.57 ) will be, we set £ = and get 



ndt = 0.8 x 10~ n s~ 2 |<5m 2 v | 5 / 6 c 17//6 , which means that the change in 5P Z 
between t res and t\ is negligible. So we can immediately give up the static 
approximations for (p~ and 



4.7.2 Mathematical Pendulum 

After the quasi-static approximation for (p~ breaks down, the first term in 



(|4.28| ) will be the dominant one due to the rapid growth of 5P Z , so we neglect 



the other terms. In ( |4.24j ), we can neglect the terms with the slowly varying 
Since b A 3> b s -& + sin<£> + , P~ will have a much larger effect on (p~ than on 
d(5P z )/dt, so that soon 3> P~, and we can thus neglect the third term 
in Q4.24j) . Therefore we get the simplified equations 

-jbPz ~ b s (fi + siny2 + ) sinyj" (4.58) 
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j t ip~ « b A 6P z . (4.59) 

Together, they give the second order equation 

(p~ ps (7 sin <y9~, (4.60) 

where (7 = bAb s ('d + siny) + ). We have now found a very simple description of 
the system. In the next two subsections, we will discuss the features of this 
equation. Afterwards, we discuss the consequences for our system. 

We take t — to to be the initial time. Then the initial conditions are given 

by 

= ^ ~ b ^P s ) b < 1, yT^O. (4.61) 

For defmiteness, we take cp^ to be positive, which implies 5P Z > 0. 

Let us first assume that g is constant. Then the differential equation cor- 
responds to a mathematical pendulum, where g is the acceleration of gravity, 
see Fig. |4l^. Usually, we could apply the small angle expansion around the 
stable minimum of the potential energy. However, in our case ip~ = corre- 
sponds to the meta-stable maximum of the potential energy. So our system 
will perform large-amplitude oscillations. We denote the amplitude of the os- 
cillations with <^max; so small '/'max corresponds to large amplitude. For illus- 
tration, we can define the analogy of a potential energy E pot = — g(l — cos<^~) 
and a kinetic energy E^ Q = |(<^~) 2 . Conservation of energy implies that the 
total energy E tot = E pot (t ) + E kin (t ) = const. Since E tot (t ) < 0, we know 
that the system will oscillate around the stable point (p~ = n with constant 
amplitude y^max = \(po \, and that ip~ will never pass the meta-stable point 
ip~ = 0. We can also state that due to the non-linearity of sin (p~ , the oscilla- 
tion frequency of the system, u, will be smaller than the oscillation frequency 
yRj/2-K for the linear small-angle approximation. In fact, for E tot — > 0, the 



oscillation frequency goes to zero. In Fig. [4.4| , we have plotted the ratio 
2itv/ Jg~ as a function of the amplitude v?max- 

Another important fact is that the time average of cosy? - , (cosy?~), is 
greater than for small y^max^ i.e. for large amplitudes. This results from 
the fact that the system develops relatively slowly close to the turning points 
±y?~ ax , where the kinetic energy is small. We have plotted (cosy?") as a 



function of the amplitude y?~ ax in Fig. |4T5 . 

Next, we discuss the differential equation with time dependent g. It is 
sufficient if we deduce the effect on the amplitude <^m ax ; we can then use v?max 



and Figs. [O] and [15] to derive the oscillation frequency and (cosy? ). 
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la 




Figure 4.3: Pictorial view of our pendulum. 




iogiof^max) 



Figure 4.4: Frequency of the pendulum as a func- 
tion of the amplitude. 
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iogloOmax) 



Figure 4.5: The average cos ip as a function of the 
amplitude. 



We will assume that g is approximately constant over the time of half 
the oscillation period, r/2, where r = 1/v > 2n / \fg. Then it is easy to 
calculate the change in the amplitude dip~ during the time dt = r/2. Taking 
the differential of the potential energy, we get 



-dg (1 — cos<^~ 



g d(cosip 



(4.62) 



Since the acceleration of gravity, g, is no longer constant, the total energy is 
not conserved, so that 



E tot = E kin + E pot = ip (p -#siny2 (p - g{l - cos</? ) 

= -g{l-cosif-), (4.63) 

where the first two terms cancel due to the differential equation ( |4.60| ). We 
get the differential 

dE tot w -dg (l - (cos^ - )) , (4.64) 
where we have time averaged cos ip~ . Equating ( 4.62 ) and ( [4.64 ), we get 



d(cos ip~ 



dg ( 



(4.65) 
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If we use dg = g dt = g r/2 and d(cos \ip |) = d(cos \ip | ) = — sin 1 99 \d\ip \, 
we finally get 

%maxl = ? T / 2 _ 1 ( C OS<^ max - (cOS^ - )) (4.66) 



fl , sm|^ maxl 

as the change in the amplitude during one half oscillation, where of course 
cos<^~ ax > (cosy? - ). We see that g > will result in a decreasing ampli- 
tude, i.e. increasing |v? max |. Analogously, g < will result in an increasing 
amplitude. If thereby |y? max | + ^|y? max | < 0, E tot > 0, so the system passes 
the meta-stable point and accelerates on the other side instead of turning. 
Then it is better to use ( [4.64[ ) instead of the differential for the amplitude to 



describe how the system evolves. 

Our approximation breaks down if |v? max | + ^|y? max | ^ 1> since then 
(cosf max — (cos<^ - )) — > 0, sin -1 <^ max — > oo and r — > oo. Also, g = const 
might not apply any longer. 

Now we can already make a statement about the neutrino asymmetry, 
represented by SP Z = (p~ /h^. If E tot < 0, the asymmetry will change sign 
at the turning points ±y? max , i.e. the asymmetry oscillates. If E tot > 0, the 



asymmetry will oscillate between yj2E tot /bA and J2(E iot + g)/bA, but will 
not change sign. 

4.7.3 The Factor g 

The next step is to derive g. The factors 6a and b s in g only vary slowly with 
time, and we will take them to be constant. Thus the main contribution 
to g comes from $ + sin ip + . For convenience, we define two new variables, 
a = $ + sin ip + and (3 = $ + cos (p + . Then g = 5 s &a«, so that g oc a. We get 

a = -Da + (b c -b T )(3 (4.67) 
ft = -Df3- {b c - b T )a + b s cos (4.68) 

where we have neglected all terms dependent on $~ . 

Initially, (p~ <C 1, so we will first consider the simpler case where cos<£>~ = 
1, sin<£> - = 0. Then we can apply the quasi-static approximation for $ + and 
ip + . Thus, 

b s {b c -b T ) Q . 
°* = D 2 + (be - brf (469) 
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Now we can determine g = b s b^a qs with the help of 

„ b s (-b T Ab 2 + D 2 Ab + D 2 b c ) 

= H [Di+Avr ' (471) 

where Ab = b c — fey. 

For small Ab, a qs is positive. However, a changes sign when the numer- 
ator becomes zero, i.e. 



Using b c ~ 6-r = D/k D , we find that the term bcbx/D 2 A;^ 2 ^> 1. Therefore, 

Ab a ^ D = k D b T . (4.73) 
This result can be expressed in terms of the dimensionless variable £ = 



1 - T(t)/T. Then we get f a « fc D /6, see (|4T55|) . 

In summary, g will increase at first and then decrease for Ab > Ab a . We 
have plotted the parameter-independent dimensionless variables a qs c/s and 



^a qs c/s = d qs T res c/ HTs as a function of £ in Figs. |16l and |4~7j, respectively. 
We can see that d qs is constant to lowest order for very small £. Therefore, 
we expand g in terms of £, 



flf = 6 



i + 0(e)- (4-74) 

J- res 



This approximation is valid for £ < £ Q . 

Now we consider the situation where changes in time. We assume 
that the oscillation period r = 0(1/ y/g) is much smaller than the damping 
time scales 1/D and 1/Ab. Then we can approximate sin<£>~ and cos<£> - by 
their time average values, and (cosy? - ), respectively, when describing the 
evolution of a and (3. Thus 

a w -Da + Abp (4.75) 
/? « -D/3- A6« + 6 s (cos^ - ). (4.76) 



Note that here we can easily neglect i? , since i? cx (sin </> ) = 0. The coeffi- 
cient Ab in (|4.75|) and Q4.76|) induces an oscillation of the two variables, which 
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0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 



Figure 4.6: Evolution of -a qs , given in ( |4.69| ), as 
a function of £ = 1 — T/T res . 

a qs [10+ 3 ] 




0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 



Figure 4.7: Evolution of f^a qs = f^f^qs as a 
function of £. 



43 



will be damped by the coefficient D toward their quasi-static approximations, 
given by 



a* = wtmA co ^)- ( 4 - 78 ) 

4.7.4 Proof for Oscillations 

It is now easy to describe the evolution of our system. We will restrict 
ourselves to proving that the pendulum, and therefore 5P Z , will oscillate for 
a certain range of mixing parameters. 5P Z will change sign at least once if g 
is positive during the first oscillation of <p~ , since then the amplitude ^"ax- 
decreases. 

We start with the initial values given in ( |4.61| ). For small tp~ 



linearize the differential equation (|4.60D , so that 



(p~ « g(p~. (4.79) 

Furthermore, we linearize g as we have done in (|4.74|) and substitute dt by 
as shown in ( [4.561) . Then we get 



= h&~, (4.80) 

where h = 1.02 x 10 2 S 2 c~ 4 / 3 |£m 2 v | 2 / 3 . The solutions for this equation are 
the Airy-functions. 

When (p~ = 0(1), the linear approximation breaks down, and (p~ will 
start oscillating with a frequency u <^ ^/g/2n. We want g oc a to be positive, 
so this first oscillation has to start before £ a . Thus we demand < £ Q , 
where is defined by (p~(£<p-) = 1. We find numerically that we can fit 
the solution of this condition with 

s 3A \Sm 2 cV \ > 10- 15 - 1 (4.81) 

for \Sm 2 \ = 10" 9 -10 5 eV 2 with an error smaller than 3%. 

Of course, we need to take into account that the quasi-static solution for 
a given in |4.69| breaks down, since (tp~) < 1 during oscillations. However, as 



we can see in ( |4.75|) and (|4.76|) , a does not depend on (cos^ ) directly, but 
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indirectly through Ab(3. Thus, the change in (cos</>~) will first affect at after 
a time scale 1/Ab. So if the condition Ab <C ^fg holds at £ = then a 
stays positive during the first few oscillation periods, and thus the amplitude 
will decrease. 

Since Ab/ \fg oc y/£, we can use the stronger condition Ab/^/g <C 1 for 
£ = £ Q . Then we get 

s|5m 2 v |- 1/6 > 3 x 1(T 6 . (4.82) 

We have now found the conditions for which 5P Z oscillates. 

We can also estimate the amplitude with which the neutrino asymmetry 
SP Z oscillates. An upper bound is given if we use the maximal value for g, 
i.e. at £ = £q,, and assume E tot = 0. Then 5P Z = (p~ jb& is maximal at 
<p~ = Ti, i.e. E kin = —E pot = 2g. We derive 

(<5-Pz)max 




1.4 x 10-*s\6miv\ 1/b c- 6 ' b . (4.83) 



We can compare this result with the numerical solution given in Section 
[4.2| . Numerically we find during the period where SP Z oscillates that 5P Z < 
7.53 x 10~ 7 for the parameters Sm 2 = —1 eV 2 , s = 10 -4 . This is in very good 
agreement with our analytical value (SP z ) maiX = 1.4 x 10~ 6 . The difference is 
due to the fact that when g is maximal, the amplitude of tp~ is not maximal, 



and thus a factor of \pl has to be replaced by y 1 + cosy? max . Furthermore, 
a is slightly decreased due to (cos<£> - ) < 1. 
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Chapter 5 
Conclusions 



We have analytically examined neutrino oscillations between a sterile and 
an active neutrino in the Early Universe, neglecting the neutrino momentum 
distribution. Our main achievement has been to prove that this system, 
which is known to create a large neutrino asymmetry, exhibits oscillations of 
this neutrino asymmetry for a large range of mixing parameters, as shown 
in Fig. |5.1| . We conclude that these asymmetry oscillations, which have been 
encountered in only some of the numerical calculations, can not arise from 
numerical instabilities. 

With the methods presented in our work, the complete analytical treat- 
ment of this system seems to have become feasible. Naturally, the next step 
will be to include the effects from the neutrino momentum distribution on 
the oscillatory behavior. Then the analytical approach will clearly be supe- 



rior to the numerical one, which continues to encounter many difficulties [49 



Finally, the duration of the asymmetry oscillations and the power law of the 
asymmetry growth after the oscillations have ceased should be derived. 

Understanding the mechanism of u a -i/ s oscillations in the Early Universe 
is very important. On the one hand, if the primordial abundances are deter- 
mined with sufficient precision, it becomes possible to constrain the mixing 
parameters of the v a — u s system. As a result, some of the models which 
have been proposed to explain the current experimental situation could be 
excluded. 

On the other hand, if future neutrino experiments, such as MiniBooNE 



42| , prove the existence of a sterile neutrino, it is necessary to understand 
the impact of the v s on the primordial abundances to test whether BBN is a 
consistent theory. 
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log 10 (\5mly\) 




log lo [sin^20] 



Figure 5.1: Bounds on the parameter space derived in our work. In the 
central white region, our proof of asymmetry oscillations is valid. In the dark- 
shaded region, the sterile neutrinos come into equilibrium, see ([4.51 ). In the 
light-shaded region, the factor g starts decreasing again before ip~ = 0(1), 
see ( 1.81 ) ; we cannot say whether the asymmetry oscillates in this region 



of parameter space. Below the lower left line, 5P Z does not reach its quasi- 
static solution before resonance, see ( |4.47|) . Our small-angle approximation, 
s <C kn, is valid only to the left of the thick dashed line. To the left of 
the thin dashed line, the condition ( f4.82|) does not apply; still, there might 
be oscillations. Above the horizontal thick line, the universe would be over- 
closed [g. 
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Appendix A 

Quasi-Static Approximation 
and its Range of Validity 



The quasi-static approximation of a variable x which follows the differential 
equation 

^x(t) = -D(t)x(t) + f(t) (A.f ) 

is given by 

Xqs = Jj- (A- 2) 

If / and D depend on time, x qs also changes its value. Therefore, x qs will 
only be a good approximation if the quasi-static condition is fulfilled, i.e. 

\ x qs\ < Dx qs . (A. 3) 

We have derived all values r x = x qs /x qs for our system of differential equa- 
tions Q4.2ij| )-(4.28) assuming that all variables except for P+ are at their 
static values, where H = —T/T. 

rtt = _ = _ 6i/ + (A.4) 

r,+ - Jkr—BH-^ . ° (A.5) 
sin (p+ b c - b T D 2 + Jp c _ bx y 



Here we have divided by sin<^+ instead of since the differential equation 
for ip + has the form jjix(t) = —D(t) sin[x(t)] + f(t). In the case of cp~, we 
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can neglect the sin as long as (f < 1. Furthermore, 



(^-)qs el +6* 



K^)qs e 1 



(5P 2 ) qs e 1 + e' 



qs 



-H -8-6- — V- 1-3 



& c A (6 C -6r) 2 



6 C - 6t V ^ 2 + ( & c - &r) J 



?V = — = -6ff - &c + r p - 



where 



D 2 + (6 C - 6 r ) 2 

e = 



2P+b A (b c -b T y 

i = - g f 2 + 6-Ml-2 (&C ~' T)2 
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